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^ . Abstract 

We calculate the electric screening mass in hot hadronic matter using two dif- 
ferent approaches, chiral perturbation theory and the relativistic virial expansion 
with empirical phase shifts, and compare the results to each other and to a gas of 
free pions and p mesons. We also compute the electric screening mass for noninter- 
acting, charged bosons with mass m on a lattice to study likely finite size effects in 
lattice gauge theory simulations of continuum QCD. For a lattice of given size, the 
continuum can be properly represented only for a window in the ratio T/m. 
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1 Introduction 



The advent of a new generation of heavy ion accelerators is making it possible to 
produce hadronic matter at high temperatures and densities. Specifically, we refer 
to temperatures near the pion mass and densities greater than the density at the 
center of atomic nuclei. Little is known about the properties of these hot and dense 
systems and of their expected transition to a phase consisting of quarks and gluons. 
In principle, both of these phases are described by QCD, but the structure of this 
theory is especially complicated for the above mentioned region of temperatures and 
densities. It is therefore reasonable to use a wide array of techniques to investigate 
different aspects of the behavior of hot and dense hadronic matter. 

One way in which we can investigate the properties of this many-body system 
is to study its response to a small perturbation using linear response theory. This 
response can then be expressed in terms of correlation functions unperturbed by the 
presence of the probe. For example, we could ask for the response of an electrically 
neutral system to an applied, static, electric field. The system responds to a weak 
perturbation, such as a heavy, charged lepton or hadron, by dynamically Debye 
screening the long range Coulomb force. The Debye screening length is independent 
of the external perturbation. 

In Ref. [P the inverse of the Debye screening length, the electric screening 
mass, niei, was studied in hot QCD. Interestingly, this can be done in two different 
ways. One way is to compute the photon self-energy. The electric screening mass 
squared is then the static, infrared limit of the time-time component of the self- 
energy. Another way is through an identity which relates the electric screening mass 
squared to the pressure 



where /i is the electric charge chemical potential. Both ways of computing the 
screening mass are exact. 




(1) 
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In this work, we will use Eq. (|^) to calculate the electric screening mass for a hot 
pion gas. This is done using three very different techniques, namely, the relativistic 
virial expansion, chiral perturbation theory, and lattice theory, which are discussed 
in successive sections. In the virial approach, dynamical information obtained from 
empirical two-body scattering phase shifts is used to compute the pressure for an 
interacting pion gas with a non-zero chemical potential and hence, from Eq. (|^), 
the electric screening mass. In the following section the screening mass is computed 
using finite temperature chiral perturbation theory, extended to the case of a non- 
zero chemical potential. The results of chiral perturbation theory are found to agree 
exactly with those of the virial expansion in the low temperature Boltzmann limit, 
T « ni-,^, to order (T/m,r)^'^^- Finally, we compute the screening mass for free, 
massive, bosons on a lattice. The ratio of the screening mass on the lattice to that 
in the continuum is studied parametrically as a function of T/m, where m is the 
mass of the boson. This gives an indication of how large a lattice is needed in order 
that lattice gauge theory properly approach the continuum QCD limit. See Fig. 1, 
taken from Ref. Q. 

Each of the above mentioned methods, as might be expected, has its advantages 
and disadvantages. The relativistic virial expansion demonstrates how the influence 
of both resonant and repulsive interactions may be included in calculations of the 
electric screening mass at temperatures close to the pion mass. Since the virial 
expansion can be expressed as an expansion in powers of the density, the results of 
this approach very likely contain the right physics for dilute systems. If, however, 
the system is dense, three- and higher-body interactions are significant. Extracting 
this information from the empirical phase shifts is difficult. The chiral perturbation 
theory approach is of interest since it contains many of the low energy properties of 
QCD. It also explains some of the low energy hadron phenomenology successfully. 
Furthermore, since the Lagrangian is known, many quantities of physical interest 
can be studied. A limitation of this approach, though, is that resonant interactions 
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are not fully accounted for. These may be expected to contribute significantly for 
temperatures close to the pion mass. Currently, lattice gauge theory is a popular 
technique to understand the structure of strongly interacting matter. This includes 
studies of various correlation functions of mesons on the lattice. While lattice 
gauge theory is in principle very powerful, finite size effects are important. Analytic 
calculations for free, massive, bosons on the lattice are therefore very useful in 
quantifying the sizes of these effects. 

Besides being of intrinsic physical interest, our calculation of the electric screen- 
ing mass is also illustrative because these techniques may be used to compute dis- 
persion relations and other response functions. The importance of alternative tech- 
niques to calculate various correlation functions has been discussed by Shuryak 
who has used both experimental phase shifts as well as QCD sum rule techniques 
to compute dispersion relations for hot hadronic matter Q. In this way we may 
get a better physical understanding of the immense information that may in prin- 
ciple be available from both heavy ion experiments and future lattice gauge theory 
simulations, especially with the Teraflop project. 

2 Relativistic Virial Expansion 

The relativistic virial expansion has been used recently to compute the thermo- 
dynamic properties of a dilute gas of interacting hadrons @, This approach is 
appealing because it allows one to systematically include the effects of both reso- 
nant and repulsive interactions in finite temperature hadronic matter by relating 
the state variables to the known, empirical, phase shifts. In this section we will use 
the relativistic virial expansion to compute the electric screening mass for a gas of 
interacting pions. It will be shown in the following section that the virial expansion 
results also provide an excellent check of chiral perturbation theory calculations 
with nonzero chemical potentials. 
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The relativistic virial expansion, introduced by Dashen and co-workers Q, re- 
lates thermodynamic state variables of a system of interacting particles to the 5- 
matrix. The partition function is separable into a product of the non-interacting 
partition function Zq and an interacting term which is proportional to bi-linear 
products of the S'-matrix and its inverse. The logarithm of the partition function 
can be written as 

\nZ = \nZo + Y.z\H\^b{hM), (2) 

where Zj = ex.p{f3fj,j) for j = 1,2 are the fugacities. The virial coefficients 6(ii,i2) 
are written as 

X ■a.,j,lAS-\E)^Sm,- (3) 

In the above, the inverse temperature is denoted by /? while V, P and E stand for 
the volume, the total center of mass momentum and energy, respectively. The labels 
ii and «2 refer to a channel of the 5-matrix which has an initial state containing 
h +^2 particles-the trace is therefore over all combinations of particle number. The 
symbol A denotes the symmetrization (anti-symmetrization) operator for a system 
of bosons (fermions) while the expression with the double-sided arrow is defined as 

S-^ (d /be) S = S-\dS/dE) - {dS-^/dE)S . (4) 

The subscript c refers to the trace over all connected diagrams. The lowest virial 
coefficient B2 = b{ii,i2)/V as V ^ 00 corresponds to the case where h = 12 = 1- 

At temperatures close to the pion mass, the system is assumed to be sufficiently 
dilute for the hadrons to interact chiefly via elastic collisions. This assumption 
greatly simplifies Eq. (^) since the 5-matrix can be expressed in terms of the phase 
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shifts 5/ as 

S{E) = Y,{21 + 1)(2/ + 1) exp(25f ) , (5) 

1,1 

where / and / denote the angular momentum and isospin, respectively. The above 
assumption also implies that two-body interactions are dominant relative to the 
three-body and higher terms. The two-body interactions are expressed via the 
second virial coefficient 

The factor gi j = (2/-|-l)(2/ + l) is the degeneracy of the (/,/) channel and M is the 
invariant mass of the interacting pair at threshold. The prime over the summation 
sign denotes that for given I the sum over / is restricted to values consistent with 
statistics. If the phase shifts 5/^ — > as — > M, an integration by parts yields 

B2 = ^ dE E'Ki {(3E)Y^ gij6{ . (7) 

For further details on the behaviour of the virial coefficients and expressions for the 
thermodynamic state variables, we refer the reader to Ref. 

In the pion gas, the vr"*", tt~ and vr'^ have chemical potentials +fi, — /i and zero, 
respectively. The chemical potential in Eq. 1^, fii + fi2 = fJ-Q, corresponds to the 
net conserved charge Q in each scattering channel contributing to Eq. (|2|). The vrvr 
pressure due to interactions can then be expressed as 

= E * ; Pq' = Te^''' ^2,Q , (8) 
Q 

where for —2 < Q < 2, —2^ < /xg < 2/i. The different vrvr channels contributing 
to the second virial coefficient -62,(3 each Q can then be decomposed, with the 
appropriate Clebsch factors, into the corresponding spin-isospin channels. After a 
little algebra, the interacting tttt pressure in the above equation is finally written in 
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terms of the spin-isospin phase shifts as 
T 



jint 



/ dEE'^Ki{(3E) 2cosh(2/x/3)(5^ + 2cosh(^/3)(5^ + 36l) 



6l + 36l + 6[ 



(9) 



The total pion pressure is given by the sum of the above interacting pressure and 
the ideal gas pressure 

dpp'^ r 1 1 1 



1 



UJ 



+ 



+ 



(10) 



g(a;-/.)/T _ I g(c^+M)/r _ 1 g'^/T _ j 

The electric screening mass for the interacting pion gas can then be obtained 
from Eq. (1). We get 

1 

K2WnmT,] H 

rn 

in=l 



2 e^i-nl 

miel = — 



1 POO . 

^ K2{l3nm^) + / dEE^K,{^E) {56^ + 36l 



(11) 



mlnT J 2m, 

The second term in the above equation is the contribution to the electric screening 
mass from the interactions. In Ref. Q it was shown that the pressure of an inter- 
acting pion gas was nearly identically the sum of the pressures of an ideal gas of 
pions and (Maxwellian) p mesons up to rather large values of the temperature. This 
was due to a near exact cancellation of the isospin weighted sum of the spin-zero 
phase shifts leaving only a contribution from the jj-resonant phase shift. Unlike 
the pressure, however, the screening mass for an interacting vrvr gas is not the sum 
of the screening masses of an ideal gas of vr and p mesons. In Fig. 2 we plot the 
electric screening mass squared (in units of e^T^) for an ideal gas of vr's and /j's as 
well as for an interacting pion gas using the expression in Eq. (11). They agree fairly 
well at lower temperatures but disagree by 10% or more at temperatures above the 
pion mass. The screening masses in the two cases differ because of the presence of 
a repulsive 5q piece in the interacting vrvr screening mass in Eq. (11). In its absence, 
they would agree almost exactly for the temperatures considered. 

In the above formulae, the virial expansion has been truncated at the level of the 
second virial coefficient; only two-body collisions were considered. At temperatures 



above the pion mass, three- and higher-body interactions wih begin to contribute 
significantly. Since it is virtuahy hopeless to expect to obtain complete experimental 
data on the 5-matrix for m-particles in and n-particles out, the extension of this 
approach to high densities is limited. One would have to rely on models of the 
many-body interactions, as obtained from effective Lagrangians, for example. 



3 Chiral Perturbation Theory 

In this section, we will use the method of effective chiral Lagrangians ||7[ to 
calculate the electrostatic screening mass in a pion gas at finite temperature. We 
will generalize the calculation of the pressure of a pion gas performed in Ref. |^] 
to the case of a finite chemical potential associated with electric charge and then 
obtain m^j using Eq. (^. 

The partition function is given by a Euclidean functional integral 

Tr exp{-H/T) = J [dU] exp J d^xLeff^ (12) 

where U{x) = exp(i(f)"' {x)t'^ / F) is an SU{2) matrix comprising the pion field </>(x). 
The integration should be performed over all configurations which are periodic in 
Euclidean time, U (x, X4 + f3) = U (x, 2:4). The effective Lagrangian -Z^e// is expressed 
as an infinite set of terms with increasing number of derivatives or quark masses, 

Le// = L(2)+LW + L(6) + ... . (13) 

The leading term corresponds to the nonlinear u-model, 

L(2) = 1^2^ (d^U^d^U - mliU + C/t)) . (14) 

The coupling constant F is the pion decay constant in the chiral limit and tjiq is 
the pion mass in the lowest order in quark masses, 

ml = -^^("T-M + md) <uu + dd> . (15) 



Neglecting isospin breaking in the quark masses, the next order Lagrangian may be 
written in the following form Q, 

1 2 1 

= --ii[TT{d,,u^df,u)\ --i2TT[d^u^dM)TT{d^u^d,u 

+ -Uml Tr [d^Wd^u) Tr (u + 



-^{h + l4)mt[T:t{U + U^)f , (16) 

which involves four coupling constants li,...,!^. The quadratic (p terms in L^'^^ 
describe free mesons of mass niQ. The terms of higher order in in L^^^ and higher 
order Lagrangians are considered as perturbations. 

Inclusion of a finite chemical potential ^ related to the electric charge is equiv- 
alent to coupling the pions to an external constant electromagnetic vector potential 
with only temporal component A4 = — i/i being non-zero and imaginary. This 
changes time derivatives in Eqs. (0) and (|l6|) to covariant ones: 

da^" - = {das'"' - /x<5„4e°"'')/ . (17) 

This results in the corresponding frequency shifts in the imaginary time propagators 
of charged pions |jl^ , 

^± ~ 7 I ■ \2 , 2~i 2 ' (-'^^) 

[UJn ±lllY +777-0 

where u;„ = lirTn. It is important that the same shift also occurs in the vertices 
containing derivatives. 

The pressure corresponds to the nonvacuum part of the thermodynamic poten- 
tial 

p{ti) = £o-n, (19) 

n = - hm ^ln[TVe-(^-'^'3)/^], (20) 

where Eq is the zero temperature and chemical potential limit of fi. Thus, to calcu- 
late P{^) one should consider all closed-loop diagrams involving all possible cou- 
plings from Leff. We will confine ourselves here to the second order in the density 



of the pion gas and thus take into account only the diagrams with one and two ther- 
mal loops . First we consider the noninteracting pion gas thermodynamic potential 
and pion propagators and their derivatives at the origin, 



G±(0) 
DoG±(0) 



(27r)3 



(27r)3 
d^p 



n 



u;2 + {iOn ± i^Y 
u;2 + (^oj^ ± 



(21) 
(22) 
(23) 



and similar expressions for second derivatives. Here = + p^. Note that 
a single spatial derivative would give a zero result, and that DqG±{^) is zero at 
^ = 0. These terms with temporal covariant derivatives contribute new terms to 
P(/i) which are absent at = 0. The corresponding temperature dependent finite 
parts of the above expressions are 

d'^p 



(0) 



In 1 



-((^+m)/t 



)+ln( 



1 



d^p 1 



DoG±(0) 



2 O/i 



(27r)3 2u; 
S'p 



1 



g{i^+M)/T _ X 
1 



+ 



1 



(2^)2 



(Aui~p.)IT 
1 



(24) 

(25) 
(26) 



We also introduce two combinations of 5(0 (^) and 51 (/^), 



5'(/^) = 35'o(^) + 3m^fifo(/i)5i(^) , 

3 

5(^) = 35ro(/i)5'o(0) + -m^ (5(0(0)51 (^) +5'o(At)ffi(0)) . 



(27) 



To simplify comparisons with the case = 0, we follow the notation of Ref. 1^] 
in which the two- loop formula for the pressure contained functions gO) 9\ a-nd their 
combination g = 3{g'^ + mlgogi), so that 50(0) = 50, 51(0) = gi and g{0) = g{0) = g. 
In terms of goifJ-), 5i(/u), gifJ-) and g{n) the two-loop pressure at /u 7^ takes the 
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following form 



+ 



4 



967r2F4 

+ [»(») + 29(f )1 + 2(2 [9(0) + 3g{p) + 2g(l')]) 

- yI^^^(99(0) + 38<,(m)+409(m)). (28) 
Here we introduced the renormalized coupling constants P] li,...,l4 through the 



h=l^[>^+^h] (29) 



relation 

1. — I \ J , 

327r2 

where A is a logarithmically divergent term and 71 = 1/3, 72 = 2/3, 73 = —1/2 and 
74 = 2. We use dimensional regularization since it preserves gauge invariance. We 
use the mass renormalization relation 

ml = ml{l- ^[3) (30) 

so that 5^0) 515 9 and g in Eq. (p^) are functions of m^. 

It is very useful to compare the representation for the pressure P(/u) obtained 
with the method of chiral Lagrangians to the result of the virial expansion in Sect. 
2. It is of course evident that the first term in Eq. (]28|), which is the ideal gas 
pressure, is exactly the same as Eq. (|lO|). Using the low temperature expansions of 
50 and gi 

goif.) = 2cosh(^/?)^^^e-^/^(l + ^^ + ...)+0(e-2-^/^), 
g^{^,) = cosh(^/5)^;^^^e-'"-/^fl + |— + ...')+0(e-2-^/^), (31) 



(27r)3/2 V 8 m, 
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in Eq. (^) the pressure may be written as a series 

oo 

Pifi) = r(m,r/27r)3/2 J2 Bnilj) expi-nmjT) . (32) 

n=l 

The n = 1 term here is the Boltzmann hmit of the free gas pressure, Eq. ([lO|), 
while the n = 2 term contains both 0(e^'"'"/^) corrections to this Hmit and the 
contribution of the interaction, 

B2{fi) = Bi'\^^) + B',^\^,,k) . (33) 

The case of a relativistic free gas corresponds to 

BM = i?(°)(;u) = (1 + 2cosh(/./3))n-5/2 (^l + ^I- + . (34) 

On the other hand, one can start from Eq. (^), take the standard low momentum 
(g << m^r) representation for the phase shifts 

sm25{{q) = 2q^'+^ {ml + q^'^' {a{ + qH\ + ...) (35) 

and use the scattering lengths and effective radii calculated in the chiral perturbation 
theory |^] in terms of the couplings 

n 7ml ( bml - - 9 - 21, .A 



" 167r2F2 \^ 127r2F2 

A straightforward check shows that the contributions of interaction to the pressure 
agree in the two approaches and B2^^{p, li) may be written as 



2r \ ^/2 / 



vrmTr/ \ 2m7r [2 

12 



-a(/i) + m^6(/i) 



■2 



where 



a(/i) = 2 cosh(2/x/3) + 2 cosh(/i/3) Qq + + Uq , 



2 cosh(2/x/3) 6o + 2 cosh(/u/3) (6o + SaJ) + 6o + &o + . (38) 



At 



this coincides with the result obtained in Ref. 

\2 



It should be noted 

that the {dgo/d^Y terms in Eq. (28) for the pressure are proportional to sinh^(/^/3). 



They are absent at = and crucial for agreement with the virial result. 
The screening mass is finally obtained from Eq. ( p8| ) using Eq. ([l|), 



mil 



+ 



2F2 



ffi(0)5i(0) + 



247r2F4 

3 {g'm? 



/"i + 2/"2 - ^/"s + ? ) 5i(0)5i(0) 



29\ /(O) 



(39) 



16 ' V ^ ' ^ 247 5mt 

where g'l{Q) = d'^gi/d^^{^ = 0). Odd derivatives of gi with respect to /i are zero 
at /i = 0. Here we eliminated mo in favor of the physical pion mass m-,^. The first 
term is of course exactly the same as in Eq. (11). All other terms are due to the 
interaction. One can check, as in the case of pressure, that they coincide to relative 
order (T/mTr)^/^ -with what follows from the interaction part of the virial result in 
Eq. ( [Ill ) if the approximation of effective radius is used for the phase shifts. 

In numerical calculation we use the central values of the recent estimates of the 



couplings li obtained in Ref. [l^] 



h 
h 



-0.62 ±0.94, /2 = 6.28 ±0.48, 
2.9 ±2.4, /4 = 4.3 ±0.9, 



(40) 



(the coupling ^4 relates F « 87 MeV to the physical coupling = 93 MeV ) . 
The results of this calculation are presented in Fig. 3. The dashed curve represents 
noninter acting pions. The chain-dashed curve includes the contributions of inter- 
actions to order F"^ and the solid curve includes also the contributions of order 
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F"^. We have displayed these separately because the chiral perturbation theory 
is naturally expressed as an expansion in inverse powers of F. See Eqs. (28) and 
(39). There is reasonable agreement with the virial calculation up to temperatures 
of around 80 MeV. Above 100 MeV the chiral perturbation expansion does not 
seem convergent. The order F~^ result starts to blow-up. This can be traced to 
the basic derivative expansion of the Lagrangian. For example, the low momentum 
expansion of the phase shifts in Eq. (35) will be inadequate because the average 
two-pion collision energy grows with temperature. 

It should be noted that in obtaining Eq. (]28| ) and Eq. ( |39|) we have actually 
taken into account three-loop diagrams which contain up to two thermal loops. Due 
to these thermal loops the final results are quadratic in the phase space functions go 
and gi. The third, T = loop, is responsible for the renormalisation of the couplings 
li. There is, however, one three- loop diagram which does not factorize into T = 
and r 7^ loops. This is the eye-type diagram and we did not take it into account. 
As was shown in Q it gives a contribution to the pressure which is proportional to 
exp(— m7r/r), but is suppressed by an additional pre-exponential factor of T/ttIt^. It 
contributes to the 0{T'^ /m\) term in Eq. (37) and thus is related to terms in the 
expansion Eq. (^5|). This contribution should become important at T ~ m-,^ which 
explains the deviation from the virial result which uses experimental information 
on phase shifts. 



4 Free Bosons on a Lattice 

In principle Monte Carlo simulations of lattice QCD should predict the properties 
of strongly interacting matter at all temperatures, including the low temperature 
phase of hadrons. So far little has been learned about the low temperature phase 
due to limitations of finite lattice size and lattice spacing and the difficulty of doing 
calculations with light quarks. Some years ago, singlet and nonsinglet quark number 
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susceptibilities were computed on a lattice of size 8'^ x 4 [^] . These susceptibilites 
are linear combinations of the baryon and electric susceptibilities. The latter is just 
the square of the electric screening mass up to a factor of e^. Since then calculations 
have also been done on a 10^ x 6 lattice and these were used in Ref. See 
Fig. 1. The temperature range covered by this larger lattice is about 120 to 190 
MeV, which is just in the interesting regime where a crossover from hadron to quark- 
gluon degrees of freedom takes place. When comparing the results of lattice QCD 
with continuum field theory calculations it is important to have an estimate of how 
important finite lattice size and spacing effects are. This is true not only of the 
lattices just mentioned, but also for lattices to be used in the upcoming Teraflop 
project; typical lattices are expected to be 48^. Is this large enough? 

To get a handle on this question we consider a system of noninteracting, massive, 
charged scalar bosons on a lattice at finite temperature and chemical potential. 



Fermions with a chemical potential have been studied on the lattice |14| , [L5[| but 
apparently there are no reports of the boson calculation in the literature. Since this 
is a free theory, the partition function can be evaluated exactly, and from this one 
can compute the net charge (or number) density and the electric susceptibility. 

We follow here the notation of Creutz ||l^. We consider a Euclidean lattice of 
size A'^^ X Nf with equal lattice spacing a in the space and time directions. Roman 
indices run from 1 to 3, Greek indices run from 1 to 4, with 4 being the time 
direction. A lattice site is specified by x^, = ariy. The integers n have allowed 
values 

-T<"^-T- ^^^^ 

The temperature is T = 1/Nta and the physical length of a side of the 3-dimensional 
cube is L = l/NgU. Letting <I> denote the complex scalar field, the action is 



S^^4J2 ^ ^y^l jy^) +mV^cI>*(n)cI>(n), (42) 

(/,n) 
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where the first sum is over aU nearest neighbors {l,n). We shall impose periodic 
boundary conditions in the spatial directions, and of course finite temperature re- 
quires the fields to be periodic in the time direction. 

We introduce a chemical potential corresponding to the conserved charge in 
the same way as one normally introduces an electromagnetic vector potential p^ . 
That is, we make the replacement iA/^ iA4 + /i. This ensures that the chemical 
potential couples to exactly the same charge density as the time component of the 
vector potential. This maintains gauge invariance. Thus, in the action, the only 
term which changes is the following one. 



We can think of this as giving particles (going forward in time) a chemical potential 
/i and antiparticles (going backward in time) a chemical potential — /x. This is 
usually called a link hopping term. 

In order to carry out the functional integration over the fields it is convenient 
to express them in terms of their Fourier components. 



where the components of the vector k are allowed the same values as the components 
of the vector n. Inserting this Fourier decomposition into the expression for the 
action, integrating over the lattice sites with the help of 



— a 



2 [$*(n, n4 + l)e"''$(n, n^) + $*(n, n4)e-''^$(n, + 1)] 



(43) 



n 




(44) 




(45) 



we get the action in the form 




(46) 



where <I> = {(pi + i(\>'i)l\f2 and the propagator is 
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+ 



1 - - (^exp |aM + ^1 + exp |-a^ - 



• (47) 



The logarithm of the partition function is, up to an irrelevant additive constant, 
given by 



InZ = In J [d^]e-^ = -Indet 
Let us define a relativistic lattice energy e by 

i 

Let us also define the complex variable 



Then the partition function can be expressed as 

In Z = - In [2 + a^e^ - ze«'* - ^-^e""'^' 



(48) 



(49) 



(50) 



(51) 



It is expedient to differentiate with respect to the chemical potential before doing 
the summations. This gives the net particle number (or charge) of the system. 

,91nZ 



N = T- 



aTY: 



We now perform the sum over k^ analytically with the help of the formula 

HZLiM\ Nt f dz f{z) 



(52) 



K4 



c z z 



Nt 



(53) 



where C is any closed contour containing the points satisfying z^* = 1 and which 
does not include the origin z = 0. Thus 

1 r dz f{z) 



where 



2m Jc z z^t -1 



ze^-f" - z-^e-"-!^ 
2 + a^e^ — ze"-i^ — z-'^e'"-'^ 



(54) 



(55) 
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The function / has simple poles at z = exp[a(±a; — //)] where a; > is defined by 

s..h(f).|. (56) 

If one analytically continues from Euclidean space (imaginary time) to Minkowski 
space (real time), as is appropriate for obtaining a response function, then the 
Matsubara frequency 2'!TkiTi — > pQ, where po is a real, continuous energy. The 
single-particle energies are determined by the poles of the propagator. In the limit 
of a vanishing chemical potential these poles are located at pQ = ztoj where lo is 
as defined in the equation above. With some rearrangement we can replace the 
integral over the single closed contour C with integrals over three disjoint contours 
C+, C- and Co encircling the two poles of f{z) and the origin. The residues are 
easily evaluated. Dividing by the volume we obtain the number density. 



{Nsaf ^ 



n = 



(57) 



This is a familiar form. It is now easy to integrate n with respect to // to obtain the 
partition function. 

In Z = - ^ [in (l - e-('^-'^)/^) + In (l - e-('^+'^)/^)] (58) 

k 

In the continuum theory Bose-Einstein condensation occurs when the chemical 
potential approaches the mass. On the lattice the condition is slightly modified. 
The number density diverges as jU approaches the critical value determined from 

sinh (^) = ^ . (59) 

This allows for any finite number of particles on the lattice even when the tempera- 
ture goes to zero. In that limit, all particles are concentrated in the zero momentum 
mode. 

Both the partition function and the number density vanish in the zero temper- 
ature limit {Nt — ^ DO at fixed a) so long as —jJLcrit < jJi < jJ-crit- They both have 
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the correct limits as the continuum is approached: a ^ 0, Nt ^ oo, Ns ^ oo with 
Nt < Ns and T = 1/Nta fixed. 

Now we turn to the electric susceptibility. Differentiating n with respect to fi, 
and setting /i = 0, we get 

This reproduces the correct expression in the continuum limit, 



-(contmuum) = - j (,^^3 (^./^ _ ^^2 > (61) 



where here uj = y/m'^ + p^. The ratio of these two expressions, lattice/continuum, 
depends only on the single dimensionless variable T/m for a given lattice x A^£, 
since the lattice spacing can be written as a = l/NfT. 

In Fig. 4 we plot lattice/continuum for lattices of size 8^ x 4 and 10^ x 6 for 



which lattice QCD calculations have been done |12, (We recall that in the QCD 
calculations the scale is set in such a way that the p meson mass has its physical 
value. The pion is then too heavy, being about 1/2 the p mass even for the larger 
of the two lattices.) The susceptibility on the 10^ x 6 lattice gets to within 40% of 
the continuum value when T/m = 0.35. It deviates markedly for both lower and 
higher temperatures. As we shall now discuss, the deviation at low temperature is 
caused by finite lattice spacing, while the deviation at high temperature is caused 
by finite lattice volume. 

Apart from the conditions already mentioned which must be satisfied if the 
lattice is to approximate the continuum, we have another. The lattice spacing must 
be small compared to all physical length scales. Thus one must have a <C 1/m. 
This is equivalent to the condition 1/Nt <^ T/m. Therefore the departure of the 
lattice susceptibility from the continuum limit will be greater and greater as the 
temperature gets smaller. This is seen in the figure. From the expressions given 
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above we can readily evaluate the susceptibilities in the low temperature limit. 

dn 

— (lattice) ^ 2iV|^, 
— (contmuum) ^ — — j e / , 

/ J- X 7/2 

lattice/continuum = {2Txf''^Nf { — \ e™/^ . (62) 

The ratio diverges exponentially as T/m 0. To get accurate results we obviously 
cannot go too low in T/m for a fixed value of Nt- 

The high temperature limit is equivalent to letting the mass go to zero. The 
susceptibility of the lattice diverges as the mass vanishes because of the zero momen- 
tum mode. This is not true of the susceptibility in the continuum; in the continuum 
the integral is convergent in the infrared because of the factor p^dp. Therefore the 
lattice/continuum ratio also increases at large values of T/m. However, this is a 
finite lattice volume effect, not finite lattice spacing effect. We can see it in the 
following ways. If we consider a box of volume instead of the continuum limit 
then we would make the replacement 

Then the susceptibility diverges in the zero mass limit because of the p = mode, 
which was suppressed in the integral. Numerically we can see this if we increase 
the ratio Nf/Ng. To approach the zero lattice spacing, infinite volume, limit we 
require (among other conditions) that this ratio be small in order that many thermal 
wavelengths fit within the box. That is, 1/T <^ L. In this sense the lattice 8^ x 4 
is 'bigger' than the lattice 10^ x 6, as is apparent in Fig. 4. 

In Fig. 5 we plot the lattice/continuum ratio for lattices of size 48'^ by 12, 16 
and 24. These may be typical for the upcoming Teraflop project. On the whole 
the susceptibilities are much closer to the continuum values than was true for the 



20 



smaller lattices. However, the small and large T/m behavior is still apparent, as 
discussed above. The condition 1 ^ A'^t ^ A'^s is also apparent in the figure. 

Since hadrons with masses in the range of 140 to 1000 MeV and beyond are 
important for the electric screening mass in the temperature range of 50 to 170 
MeV or so, it is clear that very large lattices are required to at least reproduce the 
noninteracting gas results. Whether interaction effects, and the composite nature 
of the hadrons, are more or less sensitive to finite lattice spacing and volume is 
not known. It is known, and obvious, that any interaction effects that get big 
contributions from long wavelength modes are even more sensitive to finite lattice 
volume p^ . 

5 Conclusion 

To get insight into the nature of the expected phase transition (or rapid crossover) 
of hot hadronic matter to quark-gluon plasma, it is instructive to investigate the 
temperature dependence of different quantities characterizing the system on both 
sides of the transition point (or crossover region). We have considered here one 
such quantity, the electric screening mass, in the case of hadronic gas. We used 
two different methods: relativistic virial expansion and chiral perturbation theory. 
They give very similar results up to T ~ 80 MeV. At high temperature the appli- 
cability of both approaches become doubtful, for somewhat different reasons. In 
the virial expansion, multi-pion interactions would have to be taken into account. 
This is difficult to do theoretically, apart from the fact that multi-pion scattering 
amplitudes are practically impossible to obtain experimentally. The advantage of 
the chiral perturbation expansion is that it very generally takes into account the 
basic symmetries of QCD in an expansion in powers of the pion mass and space-time 
derivatives. This is also its limitation. The p meson cannot be generated in any 
finite order; it can arise only from an infinite series of derivative terms in the effec- 
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tive Lagrangian. As the temperature rises, the average relative energy in hadron 
cohisions also rises. Therefore, more and more derivative terms in the expansion 
must be kept. Not only is this difficult to do theoretically, but the coefficients of the 
higher order terms are not known yet from phenomenology. The application of chi- 
ral perturbation theory, at least in its present form, is also limited to temperatures 
below the pion mass. 

Numerical simulations on the lattice provide a straightforward possibility to go 
all the way from low to high temperatures. However, this approach has intrinsic 
problems due to finite lattice spacing and lattice volume effects. We have shown 
here that to reproduce the free gas results for the screening mass would require very 
large lattices. Analytic methods, such as the chiral perturbation theory and the 
virial expansion, should be useful to obtain estimates on the magnitude of the final 



volume effects for the interacting hadron gas [18|. 



^ On leave from Institute of Theoretical and Experimental Physics, Moscow 117259, 
Russia. 
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Figure Captions 

Figure 1: The square of the electic mass in units of e^T^ vs. temperature. At low 
T the two hues represent the contributions from pions, and pions plus p mesons. 
At high T the three lines represent the contributions from up and down quarks 
computed to the indicated order in the QCD coupling. The data points are from 
lattice QCD calculations on a 10^ x 6 lattice. For free massless up and down quarks 
on a lattice of this size m^j/e^T^ = 35/36 as indicated by the dashed line in the 
upper right-hand corner. In this figure the pion mass was taken to be one-half the p 
mass in order to facilitate comparison with the lattice results. Taken from Ref. 

Figure 2: Square of electric screening mass in units of e^T^ vs. temperature. The 
dashed line is the contribution from a free gas of vr mesons; the chain-dashed line 
includes free p mesons as well. The solid line is the square of the screening mass for 
an interacting pion gas computed using empirical vrvr phase shifts in the relativistic 
virial expansion. 

Figure 3: Same as Fig. 2 but from chiral perturbation theory. The dashed line 
represents free pions. The chain-dashed line includes the effects of interactions to 
order F^'^. The solid line includes the effects of interactions to order F~^. 

Figure 4: Ratio of dn/dp on the lattice of specified size to the continuum as a 
function of temperature over mass. 

Figure 5: Ratio of dn/dp on the lattice of specified size to the continuum as a 
function of temperature over mass. 
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